home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / CPPCO.z / CPPCO
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCPPPPPPPPCCCCOOOO((((3333FFFF))))                                                            CCCCPPPPPPPPCCCCOOOO((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CPPCO   - CPPCO factors a complex Hermitian positive definite matrix
  10.      stored in packed form and estimates the condition of the matrix.
  11.  
  12.      If  RCOND  is not needed, CPPFA is slightly faster.  To solve  A*X = B ,
  13.      follow CPPCO by CPPSL.  To compute  INVERSE(A)*C , follow CPPCO by CPPSL.
  14.      To compute  DETERMINANT(A) , follow CPPCO by CPPDI.  To compute
  15.      INVERSE(A) , follow CPPCO by CPPDI.
  16.  
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  19.       SUBROUTINE CPPCO(AP,N,RCOND,Z,INFO)
  20.  
  21. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  22.      On Entry
  23.  
  24.      AAAAPPPP COMPLEX (N*(N+1)/2)
  25.         the packed form of a Hermitian matrix  A .  The
  26.         columns of the upper triangle are stored sequentially
  27.         in a one-dimensional array of length  N*(N+1)/2 .
  28.         See comments below for details.
  29.  
  30.      NNNN INTEGER
  31.         the order of the matrix  A .  On Return
  32.  
  33.      AAAAPPPP an upper triangular matrix  R , stored in packed
  34.         form, so that  A = CTRANS(R)*R .
  35.         If  INFO .NE. 0 , the factorization is not complete.
  36.  
  37.      RRRRCCCCOOOONNNNDDDD REAL
  38.         an estimate of the reciprocal condition of  A .
  39.         For the system  A*X = B , relative perturbations
  40.         in  A  and  B  of size  EPSILON  may cause
  41.         relative perturbations in  X  of size  EPSILON/RCOND .
  42.         If  RCOND  is so small that the logical expression
  43.         1.0 + RCOND .EQ. 1.0
  44.         is true, then  A  may be singular to working
  45.         precision.  In particular,  RCOND  is zero  if
  46.         exact singularity is detected or the estimate
  47.         underflows.  If INFO .NE. 0 , RCOND is unchanged.
  48.  
  49.      ZZZZ COMPLEX(N)
  50.         a work vector whose contents are usually unimportant.
  51.         If  A  is singular to working precision, then  Z  is
  52.         an approximate null vector in the sense that
  53.         NORM(A*Z) = RCOND*NORM(A)*NORM(Z) .
  54.         If  INFO .NE. 0 , Z  is unchanged.
  55.  
  56.      IIIINNNNFFFFOOOO INTEGER
  57.         = 0  for normal return.
  58.         = K  signals an error condition.  The leading minor
  59.         of order  K  is not positive definite.  Packed Storage The following
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCPPPPPPPPCCCCOOOO((((3333FFFF))))                                                            CCCCPPPPPPPPCCCCOOOO((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      program segment will pack the upper triangle of a Hermitian matrix.
  75.         K = 0
  76.         DO 20 J = 1, N
  77.         DO 10 I = 1, J
  78.         K = K + 1
  79.         AP(K) = A(I,J)
  80.         10    CONTINUE
  81.         20 CONTINUE LINPACK.  This version dated 08/14/78 .  Cleve Moler,
  82.      University of New Mexico, Argonne National Lab.  Subroutines and
  83.      Functions LINPACK CPPFA BLAS CAXPY,CDOTC,CSSCAL,SCASUM Fortran
  84.      ABS,AIMAG,AMAX1,CMPLX,CONJG,REAL
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.